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1. Introduction 

Polyakov has argued in a recent series of paper that de Sitter spacetime suffers 
from an explosive production of particles, if these particles have nonvanishing self- 
interactions. He did not address quantitatively the issue of backreaction, namely, 
how the original de Sitter spacetime is changed by the produced particles. Here, 
we offer a few heuristic remarks on this issue, awaiting the definitive calculation 
of the relevant effective parameters. Incidentally, an extensive list of references on 
the problem of vacuum-energy decay can be found in the original paper J22EI anc j 
a useful follow-up paperSl (see also Refs. |H and [HI for t wo early papers). 

In four-dimensional de Sitter (dS) spacetim with a positive cosmological 
constant A, infrared higher-loop quantum effects are also expected to lead to in- 
stability of the pure (matter-less) vacuum stated With the macroscopic timescale 
being set by the inverse of the Hubble constant H^s and A being the only macro- 
scopic energy-density scale available, the produced energy density of particles is 
given by the following expression (c= h = 1): 
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where the square-root factor on the right-hand side of (flaj) has been added for 
later convenience and the dimensionless parameter 7 is a short-hand notation for 
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7v (another parameter 7m will be introduced subsequently). The parameter 7 > 
involves the scalar self-coupling constants (see further comments below) and the 
overdot on the left-hand side of (jlap stands for differentiation with respect to a 
cosmic time coordinate t to be defined later (the resulting matter perturbation 
breaks the original de Sitter symmetry). It needs to be emphasized, right from 
the start, that Polyakov's tentative result (fTa| solely relies on the careful study of 
quantum field theory in curved spacetime and does not require the introduction of 
new theories. 

The constant growth rate from (|Ta)) only holds as long as pM <C A, so that 
backreaction effects can be neglected. In order to address the further evolution 
(pm ~ A), we replace the cosmological constant A by a dynamic vacuum energy 
density pv- That is, we split the backreaction effects from the produced particles 
into a standard-matter- type component and a vacuum- type component. 

2. Dynamic Vacuum Energy 

For our purpose, the following classical actiorPD33 can b e U ged: 
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£ M = -i(90) 2 + i™ 2 2 + i«0 4 , (2b) 



where the matter part (|2b[) is entirely standard, consisting of a massive real scalar 
field 4>(x) with quartic coupling constant k > 0. In (f2a|) . ey (q) is a generic function of 
a q type field, which is essentially a nonfundamental scalar field q(x) with a constant 
nonzero expectation value qo in the equilibrium state. One particular realization of 
q(x) is by use of a three-form gauge field A(x); see Refs. RJI and fTUl for details and 
original references. 

The vacuum energy density in the corresponding gravitational field equation 
is denoted by pv and differs, in general, from the energy density ey entering the 
action. In terms of a single g-type field, the gravitating vacuum energy density takes 
the following form™. 

Pv(q)=-Pv(q) = e v (q)-q^j^. (3) 

Henceforth, we assume that pv takes sub-Planckian values and that the scalar 
particles are sufficiently heavyP 

< p v < {E P f w (2.44 x I0 18 GeV) 4 , (4a) 



m » y/pv/(E P ) 2 . (4b) 

Condition Pa|) allows us to rely on classical gravity. In the cosmological context to 
be discussed shortly, condition (|4b|) implies that the "response-time" of the particle 
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wave function is very much less than the relevant cosmological timescale. Condition 
(|4b[) in the cosmological context also suggests that the produced scalar particles are 
nonrelativistic, their rest mass being very much larger than the Gibbons-Hawking 
temperature^ of the corresponding de Sitter spacetime, m 3> Tqh = Hds/Z'x. 

The action @ is invariant under general coordinate transformations or, more 
precisely, diffcomorphisms. Diffeomorphism invariance implies energy-momentum 
conservation of the matter component (cf. Appendix E.l of Ref. 11). Here, the 
matter component is given by the dynamic vacuum energy density py — —Py and 
the scalar field 4>. 



3. FRW Equations 

Now let us restrict ourselves to a spatially flat Robertson- Walker metric (in stan- 
dard comoving coordinate d^^ with cosmic time t), a homogeneous pressureless- 
perfect-fluid component (wm = Pm/pm = 0), and a homogeneous vacuum-energy 
component (wy = Py j py = — 1). Then, the Friedmann-Robertson- Walker (FRW) 
equations [in dimensionless form given by (|6c[) and (|6dp below] are consistent only 
when using the following generalization of (|la,|) : 

pM + 3Hp M = 7 V\pv\/(Ep) 2 pv, (5a) 

and the following vacuum-energy-density equation: 

py = - 7 V\PV\/(E P ) 2 py , (5b) 

where the equation of state Py = —py has been used on the left-hand side. The 
absolute value \pv\ has been employed in the square-roots on the right-hand sides 
of ([5]), but this is, strictly speaking, not necessary if (jia|) holds. 

Before presenting the complete set of differential equations, let us comment on 
the physical interpretation of (l5al) and (|5b[) . The first set of comments concerns 
the starting point, the de Sitter-spacetime calculation giving the right-hand sides 
of (|Ta)) and (|5a| . Given the theory ©, the dimensionless parameter 7 > involves 
expressions containing positive powers of the quartic coupling constant n. The main 
results of Polyakov's investigations are, first, that higher-loop quantum corrections 
may contribute the seed for 7 (for example, from a 1-loop Feynman diagram implicit 
in the second paragraph of Sec. 5 in Ref. [2]but not evaluated explicitly) and, second, 
that a chain-reaction-type evolution can amplify microscopic effects to macroscopic 
ones [see Eq. (27) in Ref. [5] for an example of this explosive behavior and Sec. [5] 
for further discussion]. Both results are infrared effects, with the whole of de Sitter 
spacetime contributing. 

The second set of comments concerns the backreaction. Purely from the classical 
theory @, the right-hand sides of (|5a|) and (|5b[) should vanish: see, respectively, 
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Eqs. (4.6) and (4.5) in Ref . ITD1 for the case of a constant gravitational coupling pa- 
rameter G — Gn- In the stationary background of the classical de Sitter spacetime 
and with classical g-type fields entering the vacuum energy density py, the scalar 
quantum field theory corresponding to (|2b[) would give rise to particle production 
from infrared quantum effecti^l as given by ([Taj) . In the dynamic context [time- 
dependent vacuum energy density pv (t)] , this results in the expression shown on the 
right-hand side of (|5a|) . with the opposite expression on the right-hand side of (|5b[) 
from energy-momentum conservation (tracing back to the diffeomorphism invari- 
ance of the theory, including quantized-matter effects). There are further terms on 
the right-hand-sides of (I5a|) and (l5b)l . which are of the form ±jm \J\pv\l {Ep} 2 Pm 
for 7m > and which correspond to a type of stimulated emission of matter parti- 
cles (cf. Sec. 5 in Ref. d]). These terms are the outcome of the chain-reaction-type 
evolution mentioned in the previous paragraph but will, for the moment, not be 
considered explicitly (7^/ = 0). 

As promised, we give the complete set of ordinary differential equations (ODEs), 
obtained by adding the two FRW equations which follow from the action Pa|) for the 
spatially flat (k = 0) Robertson- Walker metric. Introducing dimensionless variables 
by use of appropriate powers of Ep, these ODEs are 

f v = -7 \r v \ 1/2 r v , (6a) 

r M +3hr M = + 1 \r v \ 1 / 2 r v , (6b) 

2h = -r M , (6c) 

3 h 2 — r v + r M , (6d) 

with the dimensionless Hubble parameter h[r) and the dimensionless cosmic time 
t (the overdot now stands for differentiation with respect to r). The dimensionless 
variable tm{t) > corresponds to the standard-matter energy density pjw(i) > 
with constant equation-of-state parameter wm = 0. Similarly, 7v(t) = — pv(j) 
corresponds to pv(t) = —Pv{t)- Elaborating on the previous discussion of energy- 
momentum conservation, the differential system © is seen to be consistent: ([Sc]) 
follows from taking the time derivative of (|6d[) and using (l6a|) and (|6b|) . 

In order to match the setup of (fTaj) . the boundary conditions for the ODEs (|6a[) . 
(I6bl). and (IB^l) are taken as follows: 



iy(l)=A>0, (7a) 
rAf(l)=0, (7b) 
h(l) = ^X/3, (7c) 
with the Friedmann equation (|6dj) acting as a constraint. The boundary conditions 
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Fig. 1. Numerical solution of the ODEs I0 for vacuum-energy decay parameter 7 = 1 with 
different boundary conditions Q , the upper (full) curves corresponding to an initial dimensionless 
cosmological constant A = A/(Ep) 4 = 1 and the lower (dashed) curves to A = 1/10. The curves 
are rescaled by constant numerical factors at r = 1 and by appropriate r monomials for r» 1. 
The coefficients {cy, cm, ch} used in this rescaling are defined in lllOU and take the numerical 
values {16, 12.28, 6.14}. The model universe for r ~ 1 resembles a segment of de Sitter spacetime 
ifSll with approximately constant Hubble parameter, h(r) ~ h{l) = A/3. As r increases, the 
vacuum energy density rv(r) drops from the value A to zero and the model universe approaches 
Minkowski spacetime with h = 0. 



are, therefore, characterized by the single number A = A/(Ep) 4 and the ODEs © 
by the single model-parameter 7. 



4. Exact k — FRW Solution 

For 7 = 0, the solution of the ODEs ([6]) with boundary conditions (|7|) corresponds to 
the standard (eternal) de Sitter spacetimej-^ having the following constant functions: 



(7=0) 



{Mr), r M (r), h(r)} = {A, 0, \A/3 } 



(8) 



For 7 — ¥ co, the solution rapidly approaches the standard matter-dominated FRW 
solution (see the asymptotic result below). For 7=1, the numerical solution is 
shown in Fig. [TJ 

Given the relative simplicity of the differential system (|6]) and boundary condi- 
tions |(7J), it is also possible to obtain the exact solution for 7 > 0, 



tv(t) 



r M (t) 
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h(r) 




rv(r) 



(7>0) 



(9a) 
(9b) 



h(r] 



(7>0) 



(2/7) [ (« + ey - k*] + ^1/3 [ (77 + 1) (« + 9f + (77 - 1) «"] 



V3 [ (« + e)' 7 - «•»] + (7/2) [ (77 - 1) (k + oy + (77 + 1) «>!] 



(9c) 
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with definitions 

6 = t-1, n= Vl + 12/7 2 , k=2/( 1 VX). (9d) 

Incidentally, the h solution (|9c|) results from a generalized Riccati equation, namely, 
2h + 3h 2 = rv with ry from (|9a|) . [For the case of relativistic matter (ium = 1/3) 
and the same ry from (l9al) , the relevant generalized Riccati equation is /i + 2 ft 2 = 
(2/3) ry, which can also be solved explicitly] Remark that, for fixed A > and 
t > 1, solution © reduces to © as 7 approaches from above. 

The vacuum-energy density (I9a|) is seen to drop to zero monotonically as r runs 
from 1 to infinity (the left panel of Fig. [T] shows a rescaled ry). The asymptotic 
(attractor) solution for 7 > has, moreover, an equal time-dependence for the 
vacuum and matter components: 

(4/ 7 2 )r- 2 , (10a) 

l (l + yf\ + 12/72) T- 2 , (10b) 

i (l + Vl + 12/7 2 ) r- 1 . (10c) 
The asymptotic energy-density ratio depends only on the decay parameter 7: 

lim r M /ry = - ( 7 2 + 7 vV + 12) , (11) 

which goes as (1/3) 7 2 for 7 3> 1 and as 7 for < 7 <C 1. For 7 = 0(1), the 

ratio (fTTj) is of order 1. 

5. Discussion 

Two generalizations can be mentioned at this point. First, similar results 
for the asymptotic behavior are obtained if stimulated-emission-type terms 
T7m kvl 1 / 2 r M are added to the right-hand sides of the ODEs (foa[) and (fobl) . 
but, here, we keep 7m = for simplicity. 

Second, it is also possible to consider the de Sitter perturbation ([5]) to hold for 
the spatially-closed (k = 1) Robertson-Walker metric!^ The ODEs ^) are changed 
as follows: the left-hand side of (px?)) picks up a term —2k/ a 2 and the left-hand 
side of (|6d[) a term +3 k/a 2 , where a — a(r) is the cosmic scale factor with Hubble 
parameter h = a/ a. Numerical k = 1 solutions have been obtained for the follow- 
ing boundary conditions at the "waist" of de Sitter spacetime (minimum distance 
between geodesic normals; cf. Fig. 16 in Ref. |7]): ry(0) = A > 0, rj^(O) = h(0) = 0, 
and a(0) = ^/3/A. (Sec Appendix A for corresponding analytic results.) The be- 
havior found numerically is different for 7 below or above a critical value j c which 
depends on A [for example, 7 C ~ 0.9215 for A = 1]. For 7 < j c , the k — 1 numerical 
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solution has the same asymptotics as the k — solution given by (jlOp [having, in 
particular, the same value for the energy-density ratio (jlip ]. even though different 
regions of spacetime are covered. But, for 7 > j c , the k = 1 perturbation (J5]) affects 
the de Sitter universe in a fundamentally different way: a big-crunch singularity 
is produced in a finite amount of time, just as happens for the standard matter- 
dominated k = 1 FRW universe. The underlying microscopic physics determines 
which of the two types of perturbations is more appropriate, k = or k = 1, and 
with which value of the decay constant 7. 

Let us return to the original theory ^ with 7 > 0, 7m = 0, and k = 0. 
Considering the vacuum energy density ry(r), it is then possible to estimate the 
cross-over time between the de Sitter-type behavior ((5J and the FRW-type behavior 
(|10a[) . More specifically, the "half-life" of the de Sitter universe [defined as the time 
needed to reduce an initial value ry = A from by a factor 1/2] is found to have 
the following parametric behavior from (|9ap : 



where the dimensions have been restored and the right-hand side is, as expected 
from (fTa|). proportional to the inverse of the Hubble constant (jlb[) . Remark that 
our use of the term "half-life" is not intended to imply an exponential decay, which 
is indeed not the case for (|10al) . 

The half-life (jT2"j) is a direct manifestation of backreaction effects from the 
original de Sitter-spacetime particle-production (JTa)). For 7 = 0(1), this backre- 
action timescale would be of order ~ Ep/h 1 / 2 . Such a timescale would 
be very much less than the timescale found previously)*^ which is of the order 
[(i?p) 4 /A] 1/ifas for n ~ 1 — 2 [recall the assumption that A (or py) takes a value 
far below (2?p) 4 ]. Physically, the short backreaction timescale (|T2l for 7 = 0(1) 
would be due to the chain-reaction effect mentioned in Sec. |3l This timescale of 
order l/H^s is indeed seen for the explicit solution of the simplified kinetic equa- 
tion (27) in Ref. HJ assuming the relevant overlap integrals to be of order 1. (See 
also the respective Sections 5 and 6 of Ref. [4] for other examples without and with 
chain-reaction effects.) 

Let us close with some further speculations. Assume that Polyakov's mechanism 
(jlap is relevant not only for the very early universe (T ~ Ep) but also for the 
present epoch (T ~ 3K). Then, purely phenomenologically, consider having a time- 
dependent coupling 7(7-) in the differential system If j(t) is like a step-function 
which drops to zero at a relatively recent moment (t = tf rcczc ) in the history of 
the Universe, this would fix pv(t) at later times to the constant value py(ifreeze)- 
Given that the presently observed value of pm/pv is approximately equal to 1/3, 
this would suggest that the pre-freeze value of 7 should have been of order 1. 




(12) 
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But it is also possible that the parameter 7 entering (jlap and ([5]) is significantly 
smaller than 1 and that the matter of the present Universe has a different origin 
(not solely the product of vacuum-energy decay). 

Leaving these speculations aside, the clear priority at this moment is the 
de Sitter-spacetime calculation of the effective vacuum-energy decay parameters 
7v = 7, 7m, and others if present. 

Appendix A. Series-type k — 1 FRW Solution 

The basic ODEs for perturbation ([5]) with 7>0inafc = l Robertson- Walker 
metric are 

f v + 7 \rv\ 1/2 r v = 0, (A.la) 

3 (a/a) 2 + 3 k/a 2 = r v + r M , (A.lb) 

2a/a+ (a/af + k/a 2 = r v . (A.lc) 

As discussed in Sec. [5j the boundary conditions can be taken as follows: 

rv(0) = A > , o(0) = \/3/X , 6(0) = r M (0) = . (A.2) 

With these boundary conditions, it is possible to study quantitatively the backreac- 
tion effects from vacuum-energy decay, for example, by imagining that 7 would be 
turned off for r < (corresponding to the standard de Sitter universe) and turned 
on for r > (giving the perturbed de Sitter universe). Other boundary conditions 
give similar results. 

Just as for the k = case in Sec. |H the solution of (|A.lal) with boundary 
condition from (|A.2[) can be obtained immediately, 

rv(r)=x( ■ (A.3) 

\l + h/2)(VXr)J 

With this ry(r), it remains to solve the single second-order ODE (|A.lcj) . 

Given that our interest is primarily in vacuum-energy decay for relatively small 
values of the decay constant 7 and cosmological constant A, the following Ansatz 
turns out to be useful for small enough positive values of the cosmic time r: 



a(r) = cosh(v / A/ 7 3T) 



1 + 7 E c »w (^ r )' 



71=2 



(A.4) 



where coefficient c n (^) is a polynomial in 7. Inserting Ansatz (|A.4j) into the ODE 
(|A.lc[) gives the following Taylor expansion at r = 0: 

= /o(cb) + /i(c 3 , c 2 ) r 1 + / 2 (c 4 , c 3 , c 2 ) r 2 + • • • , (A.5) 
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where the functions f n shown are linear in the leading coefficient c„+2- It is, then, 
possible to solve sequentially for the coefficients c n . The first eight coefficients are 

c 2 = 0, c 3 = -l/12, c 4 = 7/32, (A.6a) 

c 5 = (8 - 9 7 2 )/720, c 6 = 7 (-2 + 3 7 2 )/576, (A.6b) 

c 7 = - (608 - 369 7 2 + 810 7 4 ) /362880 , (A.6c) 

c 8 = 7 (1288 - 489 7 2 + 1890 7 4 ) /1935360 , (A.6d) 

c 9 = (9664- 11061 7 2 + 1296 7 4 - 17010 7 6 )/39191040. (A.6e) 

With r v (t) from (|A.3|) and o(r) from (|A.4p and (|A.6|) . the function r M (r) follows 
directly from the Friedmann Eq. (|A.lb|) . The resulting expression is also a series, 

oo 

r M (r)= 1 xY / dnh)(VXry\ (A.7) 

n=l 

where the coefficient d n for n > 2 is given by an expression involving the decay 
constant 7 and the previously determined coefficients c„, c„_i, . . ., c 2 . Specifically, 
the first eight coefficients are 

di = l, d 2 ^-3 7/4, (A.8a) 

d 3 = (-2 + 3 7 2 )/6, d 4 = 7 (6- 5 7 2 )/16, (A.8b) 

d 5 = (64- 2077 2 + 135 7 4 )/720, (A.8c) 

d 6 = -7 (28 - 377 2 + 21 7 4 )/192, (A.8d) 

d 7 = (-608 + 4317 7 2 - 3645 7 4 + 1890 7 6 ) /30240 , (A.8e) 

d s = 7 (3076 - 7770 7 2 + 4995 7 4 - 2430 7 6 ) /69120 . (A.8f ) 

The coefficients d n for 7 <C 1 and 7 > 1 are seen to effectively give alternating 
series. The same holds for coefficients c n from (|A.6[) . 

This completes the construction of a partial analytic solution of the ODEs (|A.1[) . 
leaving for the future the rigorous determination of the convergence radius of the 
series appearing in (IA.4[) and (|A.7p . Still, the qualitative behavior of the solution can 
also be obtained numerically and has already been discussed in the second paragraph 
of Sec. Comparing with the numerical solutions a num (r) and Tm, numC?") °f the 
ODEs (jA.ip for 7 < 1/10, the convergence of the series in (|A.4|) and (|A.7I) is found 
to be rather poor, with a r radius of convergence of order 1/vA. 
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Note Added in Proof 

In order to clarify the last sentence of Sec.Q] it may be helpful to replace the energy 
density py in the rest of the article by A + py and the dimensionless variable ry 
by A + ry. Then, py and pu correspond, respectively, to the vacuum-type and 
standard-matter-type energy density of the produced particles. With 7 > 0, k = 0, 
and py — pm — initially, the vacuum- type component py(t) from (|5a| drops 
monotonically to the value —A, while the standard-matter-type component pm{1) 
first increases but then drops to zero due to the expansion of the Universe. Note 
that, for the simplest possible description of the instability of de Sitter-spacetime, it 
appears necessary to have at least two types of particle-production energy densities, 
distinguished by their equation of state. 
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